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Abstract 

Wc consider linear divergence-form scalar elliptic equations and vectorial equa- 
tions for elasticity with rough {L°°{Q), O C R'^) coefficients a{x) that, in particular, 
model media with non-separated scales and high contrast in material properties. 
While the homogenization of PDEs with periodic or ergodic coefficients and well 
separated scales is now well understood, we consider here the most general case of 
arbitrary bounded coefficients. For such problems, we introduce explicit and optimal 
finite dimensional approximations of solutions that can be viewed as a theoretical 
Galerkin method with controlled error estimates, analogous to classical homoge- 
nization approximations. In particular, this approach allows one to analyze a given 
medium directly without introducing the mathematical concept of an e family of 
media as in classical homogenization. We define the flux norm as the norm of 
the potential part of the fluxes of solutions, which is equivalent to the usual H^- 
norm. We show that in the flux norm, the error associated with approximating, in a 
properly defined finite-dimensional space, the set of solutions of the aforementioned 
PDEs with rough coeflicients is equal to the error associated with approximating 
the set of solutions of the same type of PDEs with smooth coefficients in a standard 
space (e.g., piecewise polynomial). We refer to this property as the transfer prop- 
erty. A simple application of this property is the construction of finite dimensional 
approximation spaces with errors independent of the regularity and contrast of the 
coefficients and with optimal and explicit convergence rates. This transfer property 
also provides an alternative to the global harmonic change of coordinates for the 
homogenization of elliptic operators that can be extended to elasticity equations. 
The proofs of these homogenization results are based on a new class of elliptic in- 
equalities. These inequalities play the same role in our approach as the div-curl 
lemma in classical homogenization. 
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1 Introduction 



In this paper we are interested in finite dimensional approximations of solutions of scalar 
and vectorial divergence form equations with rough coefficients in Q C M.'^, d > 2. More 
precisely, in the scalar case, we consider the partial differential equation 

(-dw(^a{x)Vu{x)) = fix) x£n-fGL\n),a{x) = {aij£L^{n)} 
= on dQ, 

where i7 is a bounded subset of with a smooth boundary (e.g., C^) and a is symmetric 
and uniformly elliptic on $7. It follows that the eigenvalues of a are uniformly bounded 
from below and above by two strictly positive constants, denoted by Amin(a) and Aniax('2)- 
Precisely, for all ^ G M"^ and x G fl, 

Amin(a)|C|' < e^a(x)C < A^ax(a)|e|'. (1.2) 

In the vectorial case, we consider the equilibrium deformation of an inhomogeneous 
elastic body under a given load b G {L^{Q))'^, described by 

f-div(C(x):6(n))=6(x) xGO 
1 u = on dQ, 

where C M'^ is a bounded domain, C{x) = {Cijki{x)} is a 4th order tensor of elastic 
modulus (with the associated symmetries), u{x) G M*^ is the displacement field, and for 
ip G {Hq{Q))'^, s{ip) is the symmetric part of VV', namely. 

We assume that C is uniformly elliptic and Cijki G L°°(0). It follows that the eigenvalues 
of C are uniformly bounded from below and above by two strictly positive constants, 
denoted by Amm(C) and Amax(C)- 

The analysis of finite dimensional approximations of scalar divergence form elliptic, 
parabolic and hyberbolic equations with rough coefficients that in addition satisfy a 
Cordes-type condition in arbitrary dimensions has been performed in \b'6\ [5^ I55| . In 
these works, global harmonic coordinates are used as a coordinate transformation. We 
also refer to the work of Babuska, Caloz, and Osborn |im[H] in which a harmonic change 
of coordinates is introduced in one-dimensional and quasi-one-dimensional divergence 
form elliptic problems. 

In essence, this harmonic change of coordinates allows for the mapping of the oper- 
ator La ■= div(aV) onto the operator Lq := div{QV) where Q is symmetric positive 
and divergence-free. This latter property of Q implies that Lq can be written in both 
a divergence form and a non-divergence form operator. Using the W'^''^ regularity of 
solutions of Lqv = f (for / G L^), one is able to obtain homogenization results for the 
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operator La in the sense of finite dimensional approximations of its solution space (this 
relation with homogenization theory will be discussed in detail in section [6]). 

This harmonic change of coordinates provides the desired approximation in two- 
dimensional scalar problems, but there is no analog of such a change of coordinates 
for vectorial elasticity equations. One goal of this paper is to obtain an analogous 
homogenization approximation without relying on any coordinate change and therefore 
allowing for treatment of both scalar and vectorial problems in a unified framework. 

In section [21 we introduce a new norm, called the flux norm, defined as the L^-norm 
of the potential component of the fluxes of solutions of (jl.ip and (jl.3p . We show that 
this norm is equivalent to the usual ff^-norm. Furthermore, this new norm allows for the 
transfer of error estimates associated with a given elliptic operator div(aV) and a given 
approximation space V onto error estimates for another given elliptic operator div(a'V) 
with another approximation space V' provided that the potential part of the fluxes of 
elements of V and V span the same linear space. In this work, this transfer /mapping 
property will replace the transfer /mapping property associated with a global harmonic 
change of coordinates. 

In section [3l we show that a simple and straightforward application of the flux-norm 
transfer property is to obtain finite dimensional approximation spaces for solutions of 
(jl.ip and ()1.3p with "optimal" approximation errors independent of the regularity and 
contrast of the coefficients and the regularity of dVL. 

Another application of the transfer property of the flux norm is given in section [S] for 
controlling the approximation error associated with theoretical discontinuous Galerkin 
solutions of p.ip and (|1.3p . In this context, for elasticity equations, harmonic coordinates 
are replaced by harmonic displacements. The estimates introduced in section [5] are 
based on mapping onto divergence-free coefficients via the fiux-norm and a new class of 
inequalities introduced in section HI We believe that these inequalities are of independent 
interest for PDE theory and could be helpful in other problems. 

Connections between this work, homogenization theory and other related works will 
be discussed in section [6l 

2 The flux norm and its properties 

In this section, we will introduce the fiux-norm and describe its properties when used 
as a norm for solutions of (jl.ip (and (jl.3p ). This flux-norm is equivalent to the usual 
(r2)-norm (or (ilfg (J7))''-norm for solutions to the vectorial problem), but leads to error 
estimates that are independent of the material contrast. Furthermore, it allows for the 
transfer of error estimates associated with a given elliptic operator div(aV) and a given 
approximation space V onto error estimates for another given elliptic operator div(a'V) 
with another approximation space V provided that the potential part of the fluxes of 
elements of V and V' span the same linear space. In [53], approximation errors have been 
obtained for theoretical flnite element solutions of p.ip with arbitrarily rough coefficients 
a. These approximation errors are based on the mapping of the operator — div(aV) 
onto an non-divergence form operator —Qi.jdidj using global harmonic coordinates as 
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a change of coordinates. It is not clear how to extend this change of coordinates to 
elasticity equations, whereas the flux- norm approach has a natural extension to systems 
of equations and can be used to link error estimates on two separate operators. 

2.1 Scalar case. 

Definition 2.1. For k G (^^(0))'^, denote by kpot and k^uri the potential and divergence- 
free portions of the Weyl-Helmholtz decomposition of k. Recall that kpot and kcuri are 
orthogonal with respect to the L^-inner product, kpot is the orthogonal projection of 
k onto Lloti^) defined as the closure of the space {V/ : / G C^{^)] in {L'^{^)Y. 
kcuri is the orthogonal projection of k onto L'^urii^) defined as the closure of the space 
{e : e e (C~(17))'^ div(e) = 0} in {L^{n))'^ 

For V G H^{^), define 

llV'lU-flux := \\iayip)pot\\(L^(n)r- (2-1) 

Motivations for the flux norm 

• The (•)pot in the a-flux-norm is explained by the fact that in practice, we are inter- 
ested in fluxes (of heat, stress, oil, pollutant) entering or exiting a given domain. 
Furthermore, for a vector field ^, Jg^ ^ • nds = div(^)dx = div(^pot)dx, which 
means the flux entering or exiting is determined by the potential part of the vector 
field. Thus, as with the energy norm, ||tt||a '■= /f^(Vw)"^aVii, the flux norm has a 
natural physical interpretation. An error bound given in the flux-norm shows how 
well fluxes (of heat or stresses) are approximated. 

• While the energy norm is natural in many problems, we argue that this is no 
longer the case in the presence of high contrast. Observe that in [22, contrast 
independent error estimates are obtained by renormalizing the energy norm by 
Aniin(o). In [H], the error constants associated with the energy norm are made 
independent of the contrast by using terms that are appropriately and explicitly 
weighted by a. These modifications on the energy norm or on the error bounds 
(expressed in the energy norm) have to be introduced because, in the presence of 
high contrast in material properties, the energy norm blows up. Even in the simple 
case where a is a constant (a = aid with a > 0), the solution of (jl.ip satisfies 

/jV«)^aVn=i||VA-V||;^.(n)).- (2-2) 

Hence the energy norm squared of the solution of (jl.ip blows up like 1/a as a J, 
whereas its flux-norm is independent of a (because {aVu)pot = VA~^/) 

IKIU-flux = ||VA-V||(^2(^)^,. (2.3) 

Equation ()2.3p remains valid even when a is not a constant (this is a consequence 
of the transfer property, see Corollary 12. ip . In reservoir modeling, fluxes of oil 
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and water are the main quantities of interest to be approximated correctly. The 
energy norm is less relevant due to high contrast and has been modified (in |22] 
for instance) in order to avoid possible blow up. 

• Similar considerations of convergence in terms energies and fluxes are present in 
classical homogenization theory. Indeed, the convergence of solutions of — div(a'^Vn'^) 
/ can be expressed in terms of convergence of energies in the context of F-convergence 
[62\ [T7] (and its variational formulation) or in the terms (of weak) convergence of 
fluxes in G or //-convergence HHl EH Ea [Sg HZj (a^Vn^ ^ a°S/u°). Here, weak 
convergence of fluxes is used and no flux norm is necessary unlike in our study, 
where it arises naturally. 

Proposition 2.1. ||.||a-fiux is a norm on Hq{^). Furthermore, for all ip G Hq{^) 

Amin(a)||VV'||(i2(Q))d < ll^/'lla-flux < Amax(a) || V'i/' || (^2 (2.4) 

Proof. The proof of the left hand side of inequality ()2.4p follows by observing that 

(ViPfaVij = I {Vijf{aV7p)pot (2.5) 



J' 

Jn 



from which we deduce by Cauchy-Schwarz inequality that 

_(V^)^aVV < ||VV||L2(f,)||V^|U.fl,,. (2.6) 

□ 

The proof of the main theorem of this section will require 

Lemma 2.1. Let V be a finite dimensional linear subspace of HqIH.). For f G L^(J7), 
let u be the solution of U.l\) . Then, 

. r \\u-v\\a-anx . f WC^w - aVv)pot\\(^L^n)r .o7^ 
sup mf — —- = sup mt — — (2.7j 

Proof. Since / G L^(0), it is known that there exists w G H'^{Q) n Hq[Q) such that 

( -Aw = f X £Q , ^ 

\w = U on oil. 

We conclude by observing that for v £ V, 

\\{Vw - aVv)pot\\(L^(n))d = \\{aVu - aVv)pot\\(L^n))'i- (2-9) 

□ 
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For V, a finite dimensional linear subspace of Hq{Q,), we define 

(dwaVV) := span{div(aVi;) : v€V}. (2.10) 

Note that (divaVF) is a finite dimensional subspace of H~^{Q). 

The following theorem establishes the transfer property of the flux norm which is 
pivotal for our analysis. 

Theorem 2.1. (Transfer property of the flux norm) Let V' and V be finite- 
dimensional subspaces of //g(0). For f G L^(il) let u be the solution of lil.l]) with 
conductivity a and u' be the solution of with conductivity a'. If (divaVy) = 

(diva'Vy'); then 

sup mf = sup mf — — . (2.11) 

Remark 2.1. The usefulness of (j2.1ip can be illustrated by considering a' = I so that 
diva'V = A. Then u' G and therefore V' can be chosen as, e.g., the standard 
piecewise linear FEM space with nodal basis {fi}- The space V is then defined by its 
basis {ipi} determined by 

div(aVV'i) = Acpi (2.12) 

with Dirichlet boundary conditions (see details in section [3. l.ip . Furthermore, equation 
()2.1ip shows that the error estimate for a problem with arbitrarily rough coefficients is 
equal to the well-known error estimate for the Laplace equation. 

Remark 2.2. Equation (j2.1ip remains valid without the supremum in /. More precisely 
writing u and u' the solutions of (jl.ip with conductivities a and a' and the same right 
hand side / G L^(f]), one has 

inf ||n - 7;||a-flux = inf ||n' - v\\a'-fiux- (2.13) 
Equation (|2.13p is obtained by observing that 

\\u - v\\a-au. = IIVA-^/ + div{aVv))\\^,^^^ (2.14) 

Corollary 2.1. Let X and V be finite- dimensional subspaces of Hq{Q,). For f G L^(f]) 
let u be the solution of with conductivity a. If (divaVy) = (divVX) then 

. r \\u-v\\a-RnK . r W'^W -^v\\{L^n)r ^ 

sup mf — -—— = sup mf — — (2.15) 

Equation (|2.15p can be obtained by setting a' = / in theorem 12. II and applying lemma 

Theorem 12. II is obtained from the following proposition by noting that the right hand 
side of equation ()2.16p is the same for pairs {a,V) and {a',V') whenever div(aVy) = 
div(a'Vy)- 
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Proposition 2.2. For f G L {i}) let u be the solution of Then, 

sup mf — f-— = sup — — , (2.16) 

/GL2(Q)^;ey ||/||l2(q) ^e(divavy)^ \\^^\\{L'^{n)Y 

where 

(divaVy)^ := {z G Hl{n) : Vz; G V, {Vz,aVv) = 0}. (2.17) 
Proof. For w G H'^{Q), define 



J(u;) := ml\\{V w - aVv)pot\\(L^{p\)d. (2.18) 



Observe that 



Jiw)= inf l|Vt(^-aVu-^|Lr2m^w. (2.19) 

Additionally, observing that the space spanned by Vz for z G (divaVF)""" is the orthog- 
onal complement (in (L^(O))'^) of the space spanned by aSIv + ^, we obtain that 

J{w)= sup „f 7'^"^ . (2.20) 

Integrating by parts and applying the Cauchy-Schwarz inequality yields 

J{w) < ||Au;b2(n) sup J^ii'''^"^ . (2.21) 

ze(divaVy)^ II V2;||(L2(Q))d 

which proves 

. . \\u - t^lla-flux . Il^k^ffi) , , 

sup uif — -— < sup — — , (2.22) 

/eL2(f7)^e^ ll/llL2{n) 2G{divaVy)-L \\^H{L'^{Q.)Y 

Dividing by ||At(;||^2(Q), integrating by parts, and taking the supremum over G //^(0)n 
//q($7), we get 

J{w) (Aw,z) 
sup — — = sup sup 



106/^2 (n)ni/i(n) ll^'"^llL2(f7) ze{diva\7V)^weH'2(n)nHl^{n) l|V^II(L2(Q))d||Au;||i2(Q) 

(2.23) 

we conclude the theorem by choosing —Aw = z. □ 

The transfer property (j2.1ip for solutions can be complemented by an analogous 
property for fluxes. To this end, for a finite dimensional linear subspace V C (-L^(J7))'^ 
define 

(divaV) := {div(aC) : C^V}. (2.24) 

Observe that (divaV) is a finite dimensional subspace of H~^{Q). The proof of the 
following theorem is similar to the proof of theorem 12.11 
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Theorem 2.2. (Transfer property for fluxes) Let V' and V be finite- dimensional 
subspaces of For f G L'^(il) let u be the solution of with conductivity a 

and u' be the solution of with conductivity a! . //(divaV) = (diva'V') then 

IIKV^-C)WII(.^(.))^^ ll(a-(V^--C)WII(.^(.))^ (2.25) 

/GL2(n)Cev ll/llL2(n) /GL2(n)fe'>^' 11/11^2(0) 

Theorem 12.21 will be used in section [5] for obtaining error estimates on theoretical 
non-conforming Galerkin solutions of (jl.ip . 

Corollary 2.2. Let V be a finite- dimensional subspace of (L^(ri))'^ and X a finite- 
dimensional subspace of Hq(Q,). For f G L^(f]) let u be the solution of with 
conductivity a. If (div aV) = (divVX) then 

. JI(a(V'u-C))pot||(L2(f2))d . llVu; - V7;||(i2(f^))d 
sup mf — — = sup mf — — - — ^ — (2.26) 

/GL2(n)fev WfU^n) weH},{n)nH2{n)'"^^ \\^w\\L^{n) 

Remark 2.3. The analysis performed in this section and in the following one can be nat- 
urally extended to other types of boundary conditions (nonzero Neumann or Dirichlet). 
To support our claim, we will provide this extension in the scalar case with non-zero 
Neumann boundary conditions. We refer to subsection 18.11 for that extension. 

2.2 Vectorial case. 

For k € [L'^{Q)Y^'^ , denote by kpot the potential portion of the Weyl-Helmholtz decom- 
position of k (the orthogonal projection of k onto the closure of the space {V/ : / G 
{C^{Q)y} in (L2(rj))'^x^). Define 



Ic-flux := ||(C : e(V'))poi||(L2(n))dxd. (2.27) 
Remark 2.4. Because of the symmetries of the elasticity tensor C, one has V/ E (CQ°(r2))'^ 

(v/,(C:eWW),^,,„^,,.,= (.(/).(C:eWW)„,,„„,., (2.28) 

from which it follows that definition 12.271 would be the same if the projection was made 
on the space of symmetrized gradients. 

Proposition 2.3. ||.||c-flux is a norm on {Hq{Q))'^. Furthermore, for all G {Hq{Q))'^ 

Amin(C)||e(V')||(L2(f7))dxd < llV'llc-flux < Amax(C) ||e(V') || (^2 (j^))dxd . (2.29) 

Proof. The proof of the left hand side of inequality (|2.29|) follows by observing that 



{e{^P)f ■.C:e{iP)<\\em (L2(Q))dxd IIVjIIc-Aux- (2.30) 



m 

The fact that HV'llc-flux is a norm follows from the left hand side of inequality (|2.29|) and 
Korn's inequality [sij: i.e., for all G {H^{n))'^, 

\\^i'\\iLHn)rx'i < ^/2||e('(/')||(L2(n))dxd. (2.31) 

□ 
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For V, a finite dimensional linear subspace of {Hq{Q))'^, we define 

(divC : e{V)) := span{div(C7 : e{v)) : veV}. (2.32) 

Observe that (divC : £{V)) is a finite dimensional subspace of {H^^{fl))'^. Similarly for 
X, a finite dimensional linear subspace of (-ffg(O))'^, we define 

AX := span{At; : v G X}. (2.33) 

Theorem 2.3. Let V and V be finite- dimensional subspaces of (//g (il))*^. For b G 
(L^(J7))'^ let u be the solution of lll.3\) with elasticity C and u' be the solution of hl.3\) 
with elasticity C . If (divC : e{V)) = (divC : e{V')) then 

. f ||w-u||c-flux . f \\u' - vWc'-RuK f^„.. 

sup mf — —T^ = sup mf ——- (2.34) 

be(L2(n))d^ey ||c'||(L2(f^))d be(L2{n))di'ey' \\o\\(^L'i{p,)Y 

Corollary 2.3. Let X and V be finite- dimensional subspaces of {Hq{Q))'^. For b € 
(L^(r2))'^ let u be the solution of \L3\) with elasticity tensor C. If (divC : eiV)) = A.X 
then 



sup mf —rrrj. = sup mf — — (2.35) 



\dxd 



The proof of theorem 12.31 is analogous to the proof of theorem 12.11 
For V a finite dimensional linear subspace of (L^(r2))'^^'^ we define 

(divC : V) := span{div(C : C) : C G V}. (2.36) 

Observe that (divC : V) is a finite dimensional subspace of {H^^(Q))'^. The proof of 
the following theorem is analogous to the proof of theorem 12.11 



Theorem 2.4. Let V' and V be finite- dimensional subspaces of {L'^{Q))'^^'^ . For b G 
(L^(r2))'^ let u be the solution of liL^J]) with conductivity C and u' be the solution of liL^J]) 
with conductivity C' . If (divC : V) = (divC' : V') then 

. A\{C ■■ {e{u) - 0)pot\\(L'i(m'i-^'i . p ll(C' : (^(^) - C))pot||{L2{n))dxd 

sup mf — = sup mf — 

feG(L2(n))dCGV l|o||(L2(Q))d feG(L2(n))dCeV' l|o||{L2(n))d 

(2.37) 

Corollary 2.4. Let V be a finite- dimensional subspace of (L^(r2))'^^'^ and X a finite- 
dimensional subspace of {Il}^{Vl))'^ . For b G {L'^{Q.))'^ let u be the solution of il.3\) with 
elasticity C. If (div C : V) = (AX) then 

||(C : (e(w) - C))pot||(L2(n))'ixd . \\Vw -Vv\\,L2(n))dxd 
sup mf — = sup mf — — 

bGL2(f7)CeV l|o||{L2(n))d «>e//i(f7)nH2(f7)^ex || Au'||(^2(Q))d 

(2.38) 
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3 Application to theoretical finite element methods with 
accuracy independent of material contrast. 

In this section, we will show how, as a very simple and straightforward application, the 
flux norm can be used to construct finite dimensional approximation spaces for solutions 
of and with errors independent of the regularity and contrast of the coef- 

ficients and the regularity of 9f2 (for the basis defined in subsection I3.1.2p . A similar 
approximation problem can be found in the work of Melenk where subsets of 
such as piecewise discontinuous polynomials have been used as an approximation basis 
(for the right hand side of (jl.ip ). The main difference between [l6] and this section lies in 
the introduction of the flux- norm (||.||a-flux)) which plays a key role in our analysis, since 
the approximation error (in ||.||(j_fiux-iiorm) of the space V/j on solutions of the operator 
div(aV) is equal to the approximation error (in ||.||(i/_fiux-iiorm) of the space on solu- 
tions of the operator div(a'V) provided that div(aVV/i) = div(a'Vy^). Moreover, this 
allows us to obtain an explicit and optimal constant in the rate of convergence (theorem 
13.31 and I3.4p . To our knowledge, no explicit optimal error constant has been obtained 
for finite-dimensional approximations of the solution space of (jl.ip . This question of 
optimal approximation with respect to a linear finite dimensional space is related to the 
Kolmogorov n-width |57j . which measures how accurately a given set of functions can 
be approximated by linear spaces of dimension n in a given norm. A surprising result of 
the theory of n-widths is the non-uniqueness of the space realizing the optimal approx- 
imation [571 . A related work is also |9], in which errors in approximations to solutions 
of div(aVn) = from linear spaces generated by a finite set of boundary conditions are 
analyzed as functions of the distance to the boundary (the penetration function). 

3.1 Scalar divergence form equation 

3.1.1 Approximation with piecewise linear nodal basis functions of a regular 
tessellation of Q. 

Let Q.h be a regular tessellation of of resolution h (we refer to [12]). Let Cq be the set of 
piecewise linear functions on 0.^ with Dirichlet boundary conditions. Denote by ^pk the 
piecewise linear nodal basis elements of £q, which are localized (the support of is the 
union of simplices contiguous to the node k). Here, we will express the error estimate in 
terms of h to emphasize the analogy with classical FEM (it could be expressed in terms 
of N{h), see below if needed). 

Let be the functions associated with the piecewise linear nodal basis elements (pk 
through the equation 




(3.1) 



Define 



Vh := span{$fc} 



(3.2) 
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Theorem 3.1. For any f € L (0,), let u be the solution of Then, 

sup inf II lU-flux ^ (3.3) 

where C depends only on Q and the aspect ratios of the simplices of H.^ . 

Proof. Theorem 13.11 is a straightforward apphcation of the equation (|2.15p and the fact 
that one can approximate functions by functions from Cq in the norm with 0{h) 
accuracy (since dO. is of class C"^ solutions of the Laplace-Dirichlet operator with 
right hand sides are in i^^, we refer to [E]). □ 

Corollary 3.1. For f G L^(r2), let u he the solution of in Hq{Q) and Uh the finite 
element solution of in V^. Then, 

'^P —TTw - -\ TT^ (^-^^ 

/GL2(n) ll/llL2(n) Aminlaj 

where C depends only on and the aspect ratios of the simplices of Vth ■ 
Proof. Corollary 13.11 is a straightforward application of theorem 13.11 and inequality ([2 



□ 

Let Q be a symmetric, uniformly elliptic, divergence-free (as defined in section H]) 
matrix with entries in L°°{Q). We note that this matrix will be chosen below so that the 
solutions of divQVu = / are in if^(il) if / G L^(0) and therefore can be approximated 
by functions from £q in norm with 0{h) accuracy. It follows from [TT] that this is 
not possible for the solutions of (jl.ip . In particular, in some cases Q can be chosen to 
be the identity. 

Let be the functions associated with the piecewise linear nodal basis elements ipk 
through the equation 

f - div (a(x)V^>^(x)) = d\w{QVipk) in ^ 

Define 

:= span{cl>^}, (3.6) 

Theorem 3.2. For f e L^(f2), let u he the solution of in Hq{Q) and Uh the finite 
element solution of il.l]) in . If Q satisfies one of the inequalities of theorem \4-l\ or 
theorem \4-'A then 

\W-Uh\\Hl,in) ^ C ^ 
^^P — iTTii — TT"- ^^■^) 

feL^(n) llJllL2(n) AminlOj 

where C depends only on Q and the aspect ratios of the simplices of . 

Proof. The proof follows from the fact that if Q satisfies one of the inequalities of theorem 
l4.1l or theorem 14. 21 then solutions of — div((5Vn) = / with Dirichlet boundary conditions 
are in . The rest of the proof is similar to that of the previous corollary. □ 
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3.1.2 Approximation with eigenfunctions of the Laplace-Dirichlet operator. 

In this sub-section, we assume the minimal regularity condition (C^) on the boundary 
such that the Weyl formula holds (we refer to [l9] and references therein). 
Denote by ^'^ the eigenfunctions associated with the Laplace-Dirichlet operator in 

0, and Afc the associated eigenvalues-i.e., for € N* = {1, 2, • • • } 

" ^ (3.8) 
on dU. 

We assume that the eigenvalues are ordered-i.e., < A^+i. 

Let 9k be the functions associated with the Laplace-Dirichlet eigenfunctions ^'^ ()3.8p 
through the equation 

-div(a(x)V0fe(x)) = Afc^fc in!^ 

ek = o on an ■ ^ ' 

Here, A^ is introduced on the right hand side of (j3.9p in order to normalize 9^ 
{9k = if cl{x) = I) and can be otherwise ignored since only the span of {9k} matters. 
Define 

:=span{6'i,...,6';v(fc)}, (3-10) 

where N{h) is the integer part of \Q.\/h'^. The motivation behind our definition of is 
that its dimension corresponds to the number of degrees of freedom of piecewise linear 
functions on a regular triangulation (tessellation) of VL of resolution h. 

Theorem 3.3. For f G L^(il), let u he the solution of M.l\) . Then, 



I'lU-flux 1 / 1 



lim sup inf —777-77] = — 1=\ -j-] ■ (3.11 



f<^L^n)^^^h /i||/IU2(f7) 2V^Vr(l + 



Furthermore, the space @h leads ( asymptotically as h ^ 0) to the smallest possible 
constant in the right hand side of i3.ll]) among all subspaces of Hq{Q) with N{h), 
the integer part of\^}\/h'^, elements. 



inf sup inf%-^^^ = ^f ^11 )^(l + e(iV)) (3.12) 

v,dimiv)=N f^LHn)^ev Wfh^n) 2^Vr(l + f)iV/ ^ k v j 

where the infimum is taken with respect to all subspaces of Hq(Q,) with N elements 
and e{N) is converges to zero as N ^ 00. 

Remark 3.1. The constants in the right hand side of ()3.1ip and ()3.12p are the classical 
Kolmogorov n-width dn{A,X), understood in the "asymptotic" sense (as /i — )• for 
()3.1ip and N ^ 00 for ()3.12p ) because the Weyl formula is asymtotic. Recall that the 
n-width measures how accurately a given set of functions A C X can be approximated 
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by linear spaces En of dimension n. Writing dn{A,X) the n- width measure, it is defined 
by 

dn{A, X) := ini sup inf — (7||x 

for a normed hnear space X. In our case X = Hq{^1), A being the set of all solutions 
of (jl.ip as / spans L'^{Q) for a given a{x) and Q. It should be observed there is a slight 
difference with classical Kolmogorov n- width, indeed the flux norm ||.||a-flux used in (j3.1ip 
depends on a (as opposed to the i?Q(0)-norm ). A surprising result of the theory of n- 
widths [57] is that the space realizing the optimal approximation is not unique, therefore 
there may be subspaces, other than Qh, providing the same asymptotic constant. 

Remark 3.2. Whereas the constant in (jS.lip depends only on the dimension d, the esti- 
mate for finite h given by (jS.Sp depends explicitly on the aspect ratios of the simplices of 
0,h (the uniform bound on the ratio between the outer and inner radii of those simplices). 

Proof. Let Vh be a subspace of Hq{Q) with elements. Let {vk) be a basis of Vh- 

Let v'j^ be the functions associated with the basis elements through the equation 

( Av'i^ = -div{a{x)Vvkix)) mil. fs 13) 

\ t;^ = on dQ ' 

It follows from equation (|2.15p of theorem 12.11 that the following transfer equation holds 
sup mf — — = sup mf — — - — ^ — '—LJ— (3.14) 

where for / E L'^{Q), u is the solution of p.ip . Using the eigenfunctions of the 
Laplace-Dirichlet operator, we arrive at 

WVw-Vv'Wl^^^^). _ EZiii^^-^v'^^k? ^^^^^ 



When the supremum is taken with respect to w ^ n Hq{Q), the right hand side 
of p.lSp can be minimized by taking to be the linear span of the first [lOj/Zi*^] 
eigenfunctions of the Laplace-Dirichlet operator on 0, because with such a basis the 
first N{h) coefficients of Aid are canceled, i.e. 



with N{h) = [M/h'^]. Then 



inf sup inf ^ ^ _ 

This follows by noting that the right hand side of equation (j3.16p is less than or equal 
to ^^J^^^ and that equality is obtained for w = ^ N{h)+i- 



d 
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The optimality of the constant in translates into the optimality of the constant 
associated with Vh using the transfer equation (|3.14p . i.e. 

inf sup inf H lU-flux _ — (3.18) 

Vh,dim{Vh)=N{h) f^L2i^Q)v£Vh ||j||l2(!^) y'>'N{h)+l 

We obtain the constant in (j3.1ip by using Weyl's asymptotic formula for the eigen- 
values of the Laplace-Dirichlet operator on 0, [62] . 

,3,19, 

In equation ()3.19p . |0| is the volume of fi, d is the dimension of the physical space and T 
is the Gamma function defined by T{z) := t^~^e~* dt. It follows from equation ()3.13p 
that by defining V/j = one obtain the smallest asymptotic constant in the right hand 
side of (|3.1ip . This being said, it should be recalled that the space @h is not the unique 
space achieving this optimal constant [57] . 

For the sake of clarity, an alternate (but similar) proof is provided below. By propo- 
sition [22] 

. r \\u- -ylU-flux lkllL2(n) , . 

sup mf— iiTji = sup — — , (3.20) 

feL2(n)^^^ \\J\\l^{q) 2G{divavy)^ ^\\{L'2{n))'i 

Taking inf of both sides, we have 

■ t ■ r t^lla-fiux . r ll^llL2(n) 

mf sup mi — —r = mf sup 



V,dim{V)=Nih) f^L2(^Q-^v£V \\j\\L2{n) Vh,dim(yh)=7V(/i)2g(diyaVy)^ II V^^H (L2(Q))d 

(3.21) 

Notice that the right hand side is the inverse of Rayleigh quotient, and (divaVy)"*" 
is a co-dimension N{h) space, then by the Courant-Fischer min-max principle for the 
eigenvalues, we have 

inf sup inf ^^'"i j,,^^^""^"'' = (3.22) 

y,dim(V)=JV{h) jgi2(Q)t;GV ||/||L2(n) y/XN(h)+l 

Taking V to be Q^, then the optimal constant can be achieved asymptotically as 

□ 

Remark 3.3. Theorem 13.31 is related to Melenk's n-widths analysis for elliptic problems 
[36] where subsets of such as piecewise discontinuous polynomials have been used as 
an approximation basis. The main difference between [JB] and this section lies in the 
introduction of and the emphasis on the flux- norm (||.||a-flux) with respect to which errors 
become independent of the contrast of the coefficients and the regularity of a. Moreover, 
this allows us to obtain an explicit and optimal constant in the rate of convergence. 

Remark 3.4. Write 



\9\\H--(n) '— Z^x 

k=l 



El^(^'*'^)' (3-23) 
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Then the space also satisfies, for G [0, 1), 

hm sup inf ''^ =(^(^-r)"y^'^- (3-24) 
3.2 Vectorial elasticity equations. 

Let (ei,...,erf) be an orthonormal basis of M'^. For j G {l,...,d} and /c G N* = 
{1, 2, • • • }, let tI be the solution of 



- div (C : e{Tl)^ = CjXk^k, m 

tI = 0, on dVL, 



(3.25) 



where are the eigenfunctions (|3.8p of the scalar Laplace-Dirichlet operator in 0. Let 
M := [\Q.\/h'^] be the integer part of and be the linear space spanned by rj, 

for A; G {1, . . . ,M} and j G {1, . . .,d}. 

Remark 3.5. Eigenmodes from a vector Laplace operator work as well. We use the 
eigenfunctions for a scalar Laplace operator because they are, in principle, simpler to 
compute and because our proof uses Weyl's asymptotic formula for the eigenvalues of 
the scalar Laplace-Dirichlet operator in order to obtain the optimal constant in the right 
hand side of ()3.26p and ()3.27p . Also, the eigenfunctions for the scalar Laplace operator 
encode information about the geometry of the domain 0. 

Theorem 3.4. For b G {L?{^))'^ let u be the solution of / f73|) . Then, 



lim sup inf ,,,,,, — ^ = — ■=( ^l'*- (3.26) 

Furthermore, the space Th leads (asymptotically) to the smallest possible constant 
in the right hand side of V3.2b]) among all subspaces of Hq{Q) with 0{\Q,/h'^\) 
elements. 



inf sup inf ";;~^"^-fl- = J-f ^11 V(i + e(Ar)) (3.27) 

y,rfim(F)=7Vf,g(i2(Q))d«ey ||&||(L2(n))d 20FVr(l + f)7V/ 

where the infimum is taken with respect to all subspaces of (i7Q(Q))'^ with N ele- 
ments and e(iV) is converging towards zero as N ^ oo. 



Proof. Theorem 13.41 is a straightforward application of equation (|2.35p of theorem 12.3 
and Weyl's estimate (|3.19p (the proof is similar to the scalar case). □ 



16 



Defining ^p^ as in subsection 13.1.11 for j € {1, . . . , d} let <I>^ be the solution of 



- div {C{x) : e{<l>i)^ = CjA^k, in n, 
$( = 0, on dn, 



(3.28) 



Define 

Wh := span{^>^J, (3.29) 
Theorem 3.5. For b G (^^(O))'^ let u be the solution of i f73|) . Then, 

sup inf "l'"''"^-^"" < i^/i. (3.30) 

where K depends only on and </ie aspect ratios of the simplices ofQh- 

Proof. Theorem 13.51 is a straightforward application of equation (|2.35p and the fact that 
one can approximate functions in the norm by functions from Cq with 0{h) 
accuracy. □ 

Corollary 3.2. For b € (L^(O))'^ let u be the solution of il.3\) and Ufi the finite element 
solution of lil.3\) in Wh- Then, 

sup inf ii^ii < ^^^h. (3.31) 

where K depends only on and the aspect ratios of the simplices ofQh- 

Proof. Corollary 13.21 is a straightforward application of theorem 13.51 inequality (|2.29p 
and Korn's inequality ()2.3ip . □ 



4 A new class of inequalities 

The flux- norm (and harmonic coordinates in the scalar case [S3]) can be used to map 
a given operator div(aV) (div(C : £{u) for elasticity)) onto another operator div(a'V) 
(div(C" : e(n))). Among all elliptic operators, those with divergence- free coefhcients (as 
defined below) play a very special role in the sense that they can be written in both a 
divergence-form and a non-divergence form. We introduce a new class of inequalities for 
these operators. We show that these inequalities hold under Cordes type conditions on 
the coefficients and conjecture that they hold without these conditions. 

These inequalities will be required to hold only for divergence-free conductivities 
because, by using the flux-norm through the transfer property defined in section [2] or 
harmonic coordinates as in [53] (for the scalar case) , we can map non-divergence free con- 
ductivities onto divergence-free conductivities and hence deduce homogenization results 
on the former from inequalities on the latter. 
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4.1 Scalar case. 

Let a be the conductivity matrix associated with equation (jl.ip . In this subsection, we 
win assume that a is uniformly ehiptic, with bounded entries and divergence free-i.e., 
for all I E W^, div(aJ) = (that is each column of a is div free); alternatively, for all 

Vip.a.l = 0. (4.1) 



In 

Assume that ^2 is a bounded domain in M.'^. For a d x d matrix M, define 

d 

Hess : M := ^ didjMij. (4.2) 

We will also denote by A~^M the d x d matrix defined by 

(A-^M)ij = A-^Mij. (4.3) 

Theorem 4.1. Let a be a divergence free conductivity matrix. Then, the following 
statements are equivalent for the same constant C: 

• There exists C > such that for all u G Hq{Q,), 

(4.4) 

• There exists C > such that for all u G Hq(Q,), 

\\{diy{aV)r^An\\^,^^^<C\\n\\^,^^y (4.5) 

• Writing Oi the solutions of i TO) . For all {Ui,U2, ■ ■ ■) £ M^*, 

oo oo 

\\J2UA\\l^^.<C'Y.Ul (4.6) 



|2 

i=l i=l 



The inverse of the operator — div(aV) (with Dirichlet boundary conditions) is a 
continuous and bounded operators from onto . Moreover, foru G H~'^{Q.), 

||(divaV)-itx||^2(f,) < C||A-it/||i2(n). (4.7) 
There exists C > such that for all u G Hq{^), 

ll«lli^(o)<^^'E(dM«V-l),n) . (4.8) 
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There exists C > such that 

1 . , {Vz,aVu)L2(^n) 
7; < mi sup -r— .. .. ■ (4.9) 

There exists C > such that for all u € Hq{Q), 

II^IIl2{Q) < C||A~^Hess : {au)\\^2^Qy (4.10) 

There exists C > sitc/i i/iai /or a// u € //(^(ri), 

hllL2(f7) < <^|| Hess : (A-i(an))||^2(^^. (4.11) 



Remark 4.1. Theorem 14.11 can be related to the work of Conca and Vanninathan 
on uniform i/^-estimates in periodic homogenization, which estabhshed a similar result 
in the periodic homogenization setting. 

Proof. Let Uj £ M. Observe that 

00 00 

- div (aV J2 = (4.12) 

hence 

00 00 

- div (aV OjUj) = Y "^jUj- (4.13) 

Identifying u with Yl'jLi ^j^j^ it follows that 

. , ll^"'div(aVu)||^2(f,) 1 

mf r— — > — (4.14) 

MeL2(n) lFllL2(n) 



is equivalent to 



J2^MUn)<C'U^U. (4.15) 



Observe that equation (|4.4|) is also equivalent to 

||(divaV)-in||^,(^^ < C||A-in||i2(f,), (4.16) 

which is equivalent to the fact that the inverse of the operator — div(aV) (with Dirichlet 
boundary conditions) is a continuous and bounded operator from onto L^. Finally, 
the equivalence with (j4.8|) is a consequence of equation (|4.6p . Let us now prove the 
equivalence with equations ()4.10p and ()4.1ip . Observe that if a is a divergence free dx d 
symmetric matrix and u G H'^{^}) n Hq{^1), then 

div(aVn) = Hess : (an), (4.17) 
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since 

d d d d d 

Hess : {au) = aijdidjU + diaijdjU + djaijdiU, (4-18) 

i,j=l j=i *=i *=i i=i 

XliLi dittij = and X^^=i djaij = 0. It follows that 

div(aVn) = Hess : {au) = Hess : A~^{au), (4.19) 

which concludes the proof of the equivalence between the statements. □ 



Theorem 4.2. // a is divergence- free, then the statements of theorem \4-l\ are implied 
by the following equivalent statements with the same constant C. 

• For all u G H^{n) n H'^{n), 

l|Au||L2(n) < C\\a : Hess(n)|U2(f,). (4.20) 

• There exists C > such that for u G C^{i}), 

\\k^T{u)\\L2 <C\\k'^.T{au).k\\L2, (4.21) 
where J-{u) is the Fourier transform of u. 

Remark 4.2. Concerning equation ()4.2ip since u is compactly supported in 0, u can be 
extended by zero outside of without creating a Dirac part on its Hessian and J-{u) is 
the Fourier transform of this extension. 

Proof. Equation (j4.10p is equivalent to 

1 . , (^,A-iHess:(an))^, 

— < mf sup 77— 77— —. (4.22) 

C u^Hl{Q.) ^izL2(^Q_) \M\L'^(Q.)W\\L'^{n) 

Denoting by ip the solution of Is.-if) = if in Hl{Vl) n H'^{Vl), we obtain that (|4.22p is 
equivalent to 

1 . , (V',Hess: (a^z))^2(n) 

— < mf sup — — (4.23) 

c- u&Hi{n)^(:Hi{Q.)nm(a) \m\L-2[Q.)\\A'ip\\L-^{Q) 

Integrating by parts, we obtain that ()4.23p is equivalent to 

1 . , (a:Hess(V'),n)^2(n) 

— < mf sup 7-—— —. (4.24) 
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Since a is divergence free, a : Hess = div(aV.) and so there exists '0 such that a : 
Hess('0) = u with Dirichlet boundary conditions. For such a ip, we have 



(g : Hess(-0),M)^2^^^ ^ ||a : Hess(^-(/;;||j^2(i^) 
lkllL2(n)||A^/;||i2(n) ||A^/;||i2(f2) 

It follows that inequality ()4.24p is implied by the inequality 



(4.25) 



1 \\a : Hessl ih)\\T 2 (r)\ 

i< inf " „,,, (4.26) 

The equivalence with (|4.2ip follows from a : Hess(n) = Hess : (an) and the conservation 
of the L^-norm by the Fourier transform. □ 

Theorem 4.3. Let a be a divergence free conductivity matrix. 

• If d = 1, then the statements of theorem \4-^ are true. 

• If d = 2 and Q is convex then the statements of theorem |^.^| are true. 

• If d>2>, is convex and the following Cordes condition is satisfied 

/ (Trace[a(a;)])^ \ 
esssup,,^ [d - T^^^^^rl^^) < 1 (4-27) 

then the statements of theorem are true. 

• If d > 2, ^} is non-convex then there exists Cn > such that if the following Cordes 
condition is satisfied 

/ ( Trace[a(x)])^ \ 

esssup^gf^ I d - — f Tt \ f \] ! < (^-28 

V Trace[a-' (x)a(j;)J / 

then the statements of theorem \4-'A '^^^ true. 

Proof. In dimension one, if a is divergence free then it is a constant and the statements 
of theorem 14.21 are trivially true. Define 

/ (Trace[a(a;)])^ \ 
:= esssup,,^ [d - Trace[a^(x)a(x)] J ('-''^ 

Theorem 1.2.1 of implies that if Vl is convex and /3a < 1, then inequality (j4.20p is 
true. In dimension 2, if a is uniformly elliptic and bounded, then f3a < 1- It follows 
that if d = 2 and is convex or if d > 3, 17 is convex, and /3a < 1, then the statements 
of theorem 14.21 are true. The last statement of theorem 14.31 is a direct consequence of 
corollary 4.1 of 
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For the sake of completeness we will include the proof of three bullet points here {il. 
convex). Write C the differential operator from H'^{Q,) onto L^(r2) defined by: 

Cu ■.= ''^aijdidjU (4.30) 

Let us consider the equation 



Cu = f in Q, 
u = on (9il 



(4.31) 



The following lemma corresponds to theorem 1.2.1 of [35] (and a does not need to be 
divergence free for the validity of the following theorem). For the convenience of the 
reader, we will recall its proof in subsection 18.21 of the appendix. 

Lemma 4.1. Assume Q to be convex with C'^ -boundary. If Pa < ^ then ( |^.^j[ ) has a 
unique solution and 

II II ^ esssup^ajx) 

ImlH^nH^in) < WfUHn) (4.32) 

where a{x) := {^f=iaii{x))/ Ylij=i{aij{x) f 

Pa is a measure of the anisotropy of a. In particular, for the identity matrix one has 
Pl^ = 0. Furthermore in dimension 2 

a ■ r 2Amin(Q(a;))Amax(Q(a:)) 

Pa = l- essmf.,. + (A_(a(x)))^ ^''^'^ 

and one always have /3a < 1 provided that a is uniformly elliptic and bounded. The first 
three bullet points of theorem 14.31 follow by observing that if < 1 then 

WAmnHliu) < C|l aijdidju\\L2(^n) (4.34) 
which implies inequality ()4.20p . 

□ 



4.1.1 A brief reminder on the mapping using harmonic coordinates. 

Consider the divergence- form elliptic scalar problem (jl.ip . Let F denote the harmonic 
coordinates associated with (jl.ip -i.e.. F{x) = (Fi(a;), . . . , Fd(x)) is a d-dimensional 
vector field whose entries satisfy 

div aVFi = in , , 

(4.35) 

Fi{x) = Xi on d^. 

It is easy to show that F is a mapping from Q. onto Jl. In dimension one, F is trivially 
a homeomorphism. In dimension two, this property still holds for convex domains [Tl[5]. 
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In dimensions three and higher, F may be non-injective (even if a is smooth, we refer 

to IS], 120]). 

Define Q to be the positive symmetric d x d matrix defined by 

It is shown in [53] that Q is divergence free. Moreover, writing u the solution of (|l.ip 
and \\u\\a := Vn • aVu one has for v G Hq{^1) 

\\u - v\\a = \\u - v\\q, (4.37) 

where v := v o and u := uo solves 

-^g.,9.a,^ = ^^^oF-i (4.38) 

Note that (|4.37p allows one to transfer the error for a general conductivity matrix a to a 
special divergence- free conductivity matrix Q. Observe that the energy norm was used 
in [53] (and ()4.37p . instead of the flux norm) under bounded contrast assumptions on a. 

The approximation results obtained in |53j are based on (j4.37p and can also be 
derived by using the new class of inequalities described above for Q. 

4.2 Tensor ial case. 

Let C be the elastic stiffness matrix associated with equation (|1.3p . In this subsection, we 
will assume that C is uniformly elliptic, has bounded entries and is divergence free-i.e., 
C is such that for ah / G W^""^, div(C : /) = 0; alternatively, for all ip G {C^{9.)Y, 



(Vipf -.C -.1 = 0. (4.39) 



The inequalities given below will allow us to deduce homogenization results for arbi- 
trary elasticity tensors (not necessarily divergence-free) by using harmonic displacements 
and the flux- norm to map non-divergence free tensors onto divergence- free tensors. 

For a d X d X d tensor M, denote by Hess : M the vector 

d 

(Hess : M)k := didjUij^k- (4.40) 

Let A~^M denote the d x d x d tensor defined by 

(A-iM)ij, = A-^M,,,- fc. (4.41) 
The proof of the following theorem is almost identical to the proof of theorem 14.11 

Theorem 4.4. Let C be a divergence free elasticity tensor. The following statements 
are equivalent for the same constant 7; 
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There exists 7 > such that for all u G {HQiyi))'^ , 

M\(L^{n)Y <7||A~^div(C:e(n))||(^2(f^))d. (4.42) 

There exists 7 > such that for all u G {Hq{^)Y , 

||(div(C : e(.)))-^An||(^,(^^), < l\\u\^,^^^^a. (4.43) 

For all (C/i,C/2,...) G {W^f, 



00 d 

|2 

fc=i i=i fc=i 



where {t^} «s t/ie &asis defined in 113.25]) . 

The inverse of the operator — div(C : e(.)) (^wii/i Dirichlet boundary conditions) 
is a continuous and bounded operator from {H~'^)'^ onto {L?)'^ . Moreover, for 

u G {H-^{n)f, 

\\{divC: e(.))^^^i||(i2(f^))d < 7l|A"^n||(i2(n))d. (4.45) 
There exists 7 > such that for all u G {Hq{^)Y , 

There exists 7 > such that 

1 ((V.)^ : C : e{u))L^n) . . . 

— < mi sup -TT—r ^7— — r (4.47) 

There exists 7 > such that for all u G {Hq{Q))'^, 

\W\\iL^n)r ^ Hess : {u.C)\\^,^^y (4.48) 

There exists 7 > such that for all u G {Hq{^)Y , 

\M{LHn)Y ^ 7II Hess : {^-\u.C))\\^^^^^^^^. (4.49) 



Theorem 4.5. If C is divergence-free, the statements of theorem\4^ are implied by the 
following statement with the same constant 7. 

• For all u G {H^{n) n H'^{^)Y, 

l|An||(i2(f^))d < 7II Hess : {u.C)\\^L2^^))d- ^^■^^) 

Proof. The proof is similar to that of theorem 14.21 □ 
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4.2.1 A Cordes Condition for tensorial non-divergence form elliptic equa- 
tions 

Let us now show that the inequality in theorem l4.5l and hence the inequalities of theorem 
14.41 are satisfied if C satisfies a Cordes type condition. The proof of the following theorem 
is an adaptation of the proof of theorem 1.2.1 of [45] (note that C does not need to be 
divergence free in order for the following theorem to be valid). 

Let C denote the differential operator from {H'^{Q)'^ onto (L^(f]))'^ defined by 

{Cu)j ■.= ^CijkididkUi. (4.51) 

i,k,l 

Let us consider the equation 

{Cu = f in , , 

u = on oil. 

Let B be the dx d matrix defined by Bjm = Ylk=i ^kmkj- Let A be the dx d matrix 
defined by Aj^rn = YA,k,i=iCim.kiCij'ki- Define 

/3c := d^ - Trace[SA~^5^]. (4.53) 



Theorem 4.6. Assume is convex with a C -boundary. If Pc < 1; then |^.5^ ) has a 
unique solution and 

\M{H^nHl(n)r ^ -f'^ll/ll{L2{c))'*> (4-54) 
where K is a function of fic and ||i?^~^||(L=o(f^))dxd. 

Remark 4.3. /3c is a measure of the anisotropy of C. In particular, for the identity 
tensor, one has /Sj^ = 0. 

Proof. Let u be the solution of Cu = f with Dirichlet boundary conditions (assuming 
that it exists). Let a be a field oi d x d invertible matrices. Observe that (|4.52|) is 
equivalent to 

Au = af + Au - aCu. (4.55) 

Consider the mapping T : {H^ n H^{n) f {H^ n H^{n))'^ defined hy v = Tw, 
where v be the unique solution of the Dirichlet problem for Poisson equation 

Av = af + Aw — aCw. (4.56) 

Let us now choose a so that T is a contraction. 
Note that 

\\Twi - TwiW^H^^Himr = 11^1 - ^2||(//2nHi(n))^- (4.57) 
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Using the convexity of 17, one obtains the following classical inequality satisfied by the 
Laplace operator (see lemma 1.2.2 of |45j): 

11^1 - V2\\(H^nH^m'' < ll^(^^i - V2)\\(L^Q)y- (4.58) 

Hence, 

\\Twi - Tw2\\'^rr2nm(c,^^d <\\^iwi - W2) - aC{wi - ^2^"^ 



d d 

ej [Sjidki - ^ ajjiCij'ki)didk{w[ - W2) 

i,j,k,l=l j'=l 



2 

(4.59) 

Using the Cauchy-Schwarz inequality, we obtain that 

^ d d 

(4.60) 

( {didk{w{-w'2)?). 

i,k,l=l 

Hence, writing 

d d 

l3a,C ■= ^ {Sjlhi - ^ OijjiCijiklf, (4.61) 
i,i,k,l=l j'=l 

we obtain that 

II 1 1 2 II 1 1 2 

\\Twi - Tw2\\(^j^2^jji(^^^y < esssup^(,^l3a,c{x)\\wi - W2\\(^jj2r,Hl^(^n))d- (4-62) 
Observe that 

d d d 

f3a,c := d? - 2 ^ o-jj'Ckj'kj + ^ CY ^jj'Cij'kif- (4.63) 

j',j,k=l i,j,k,l=l j'=l 

Taking variations with respect to a, one must have, at the minimum, that for all j, m, 

d d d 

CjmkiCy^^ ajj'Cjj'ki) = y^^Ckmkj- (4.64) 

i,k,l=l j'=l k=l 

Hence, 

d d d 

5Z ^ CjmklCjj'kl = y^^Ckmkj- (4.65) 

j'=l i,k,l=l k=l 
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Let B be the matrix defined by Bjm = '^k=i^kmkj- Let A be the matrix defined by 
Aj'm = J2lk,i=i CimkiCij'ki- Then (|4.55|) can be written as 

aA = B, (4.66) 

which leads to 

a* = BA~^. (4.67) 

For such a choice, one has 

d d d 

E C£,c.)j'Cirki?= E ^U^kmkj. (4.68) 

i,j,k,l=l j'=l j,m,k=l 

Hence, at the minimum, I3a,c = l^c with 

Pc := - Trace[BA~^B'^]. (4.69) 

For that specific choice of q, if /3c < 1, then T is a contraction and we obtain the 
existence and solution of ()4.52p through the fixed point theorem. Moreover, 

ll^^^ll{L2(n))d < fh^in) + l3c\\'^u\\(^L2^Q))d, (4.70) 
which concludes the proof. 

□ 

As a direct consequence of theorem 14.51 and theorem 14.61 we obtain the following 
theorem. 

Theorem 4.7. Let C be a divergence free bounded, uniformly elliptic, fourth order 
tensor. Assume O is convex with a C'^ -boundary. If f3c, defined by |.^.53D , is strictly 
bounded from above by one, then the inequalities of theorem \4-5\ and theorem \4-4\ ^-''"^ 
satisfied. 



5 Application of the flux-norm to theoretical non-conforming 
Galerkin. 

The change of coordinates used in [53] (see also subsection I4.1.ip to obtain error esti- 
mates for finite element solutions of scalar equation (jl.ip in two-dimensions admits no 
straightforward generalization for vectorial elasticity equations. In this section, we show 
how the flux-norm can be used to obtain error estimates for theoretical discontinuous 
Galerkin solutions of (jl.ip and ()1.3p . These estimates are based on the inequalities in- 
troduced in section [3] and the control of the non-conforming error associated with the 
theoretical discontinuous Galerkin method. The control of the non-conforming error 
could be implemented by methods such as the penalization method. Its analysis is, 
however, difficult in general and will not be done here. In the scalar case, we refer to 
|52j for the control of the non-conforming error. 
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5.1 Scalar equations 

Let w G n Hq{^}) such that —Aw = f. Let u be the solution in i^Q (^7) of 

[ {V^faVu= [ {VipfVw (peH^in) (5.1) 
Jn Jn 

Let V be a finite dimensional linear subspace of {L'^{Q,))'^. 

We write Cv sua approximation of the gradient of the solution of p.ip in V obtained 
by solving (|5.2p -i.e.. is defined such that for all rj £V, 



aCv = / r/ Vw. (5.2) 

For ^ G (L^(J7))'^, denote by ^ = ^curi + ^pot the Weyl-Helmholtz decomposition of ^ 
(see Definition 12. ip . 



Definition 5.1. Write 

/Cv := sup —J— . (5.3) 

Cev IKII(L2(f7))rf 



/Cv is related to the "non-conforming error" associated with V (see for instance [19] 
chapter 10). If /Cy > then the space V must contain functions that are not exact 
gradients. Moreover, it determines the "distance" between V and L^^j (see definition 

EH). 



Definition 5.2. Write 

X>V := mi sup mr — — -7^ ^ — 

a',V' :div(a'V')=div(aV) ^'gH^ (f7)nHi (f^) II IIl^C^) 



The first minimum in ()5.4p is taken with respect to all finite dimensional linear sub- 
spaces V of {L'^{Q)Y, and all bounded uniformly elliptic matrices a' {a[j € L°°{Q)) 
such that div(a'V') = div(aV). Furthermore, u' in (|5.4p is defined as the (weak) solu- 
tion of div(a'Vn') = Aw' with Dirichlet boundary condition on dQ. Due to Theorem 
12.11 the infa'^v' :div{a'V')=div(aV) ^an be dropped. However, we keep it to emphasize the 
independence of the choice of V' and a' as long as they satisfy div(a'V') = div(aV). 

Theorem 5.1. There exists a constant C* > depending only on Amm(a) md Aniax(o) 
such that for fCy < C* , 

l|Vn - Cv||(L2(f,)). < C||/||i2(f,)(Pv +/Cv) (5.5) 



where u is the solution of h5.1\) . C,y the solution of 115. 2\} and C is a constant depending 
only on Amin(a) and Amax(a)- 
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Remark 5.1. Theorem 15.11 is in essence stating that the approximation error associated 
with V and the operator div(aV) is proportional to and /Cv is related to the non- 
conforming error associated to V. is the minimum (over a', V such that div(a'V') = 
div(aV)) approximation error associated to V' and the operator div(a'V). Hence, Py and 
the transfer property allow us to equate the accuracy of a scheme associated with V' and 
a conductivity a' to the accuracy of the scheme associated with V and the conductivity 
a provided that div(a'V') = div(aV). 

Remark 5.2. In fact, it is possible to deduce from theorem 15.11 that the maximum ap- 
proximation error associated to V and the operator div(aV) can be bounded from below 
by a multiple of (Py + ^v) (see also equation (10.1.6) of [H]). 

Remark 5.3. If the elements of V are of the form rj = ^^-^q^ l(2;eT)Vu where v belongs 
to a linear space of functions with discontinuities at the boundaries of the simplices of 
O/j then we can replace the right hand side of (|5.2|) by — J^vAw (see subsection 1.3 of 
[53]). This modification doesn't affect the validity of ()5.5p since the difference between 
the two terms remains controlled by Py. For clarity of presentation, we have used the 
formulation (|5.2p . 

In order to prove theorem 15. 11 we will need the following lemma 

Lemma 5.1. There exists C depending only on Ainin(a)5 Amax(a) such that for u G Hq{Q) 
and C G {L'^{^)Y 

(5.6) 

l|Vn-C||(L2(f^))d > ^(||(a(Vn - C))pot||(^2(Q))d + llCc«r/||{L2(n))<') (5-7) 
Proof. For the proof of ()5.7p , observe that 

- CII(L2(f^))d = ||Vn - Cpot\\{L^{n))'i + \\Ccurl\\{L^{n))''- 

Furthermore, 

II (a(Vu — C))pot II (-^2 (-Q-J^d — II {0,(\/U — Cpot))pot\ \ (/,2(-f2))d ~l" II {(^Ccurl) pot\ \ f^£^2 (^Q-j-^d 

^ '^max(o) II V'^i ''Po'll (L2(n))'' ~^ 1 1 1 1 (L^ (Q)) 



For (|5.6p . observe that 



(5.. 



{Vu - Cfa{Vu - C) = / (Vn - Cpotf{a{Vu - C))^,, + / CMVu - C) (5.9) 
It follows from Cauchy-Schwarz inequality that 

A„..n(a)||V^ - CkLHn)r ^ ^Vn - Clll^w^ " ^^^^°*"(^^(^»^ (5.10) 



+ ^max{a)\\Ccurl\\(L^{n))d 



□ 



29 



We also need the following lemma, which corresponds to lemma (10.1.1) of 

Lemma 5.2. Let H be a Hilhert space, V and be subspaces of H (Vh may not be a 
subset ofV). Assume that a(., .) is continuous bilinear form on H which is coercive on 
Vh, with respective continuity and coercivity constants C and 7. Let u &V solve 

a{u,v) = F{v) \/veV (5.11) 

where F £ H' (H' is the dual of H). Let Uh G solve 

a{uh,v) = F{v) "^vGVh (5.12) 

Then 

II II / /'I , '^A • f 1 ~ ^h,w) 

\\u — UhllH < -\ ) mi n — til H — sup 7 — 7 (5.13) 

7 "'GVh 7^6i4\{o} \H\h 



We now proceed by proving theorem 15.11 
Proof. Using lemma \57L\ we obtain that 

- CII(L2(n))<* < c(\\{a{Vu - C))pot\\^L^n))<i + \\C\\iL^n))'i'<^v) (5.14) 

Using the triangle inequality \\C\\{L2{n))'i < l|Vn - Cll(L2(f7))d + l|Vn||(i2(Q))d, we obtain 
that 

llV-u - C\\{mn))i < ^ _ fj^^ - C))pot||(i2(n))d + l|Vn||(i2(f^))d/Cvj (5.15) 

from which we deduce that 

C / 

inf \\Vu - C||(L2(n))<^ < i_c]Cy ^ev (H^^^^'' " + \Wu\\{L^{n))^^v] 

(5.16) 

(C* in the statement of the theorem is chosen so that /Cy < C* implies CfCy < 0.5). 
We obtain from lemma 15.21 that (observe that the last term in equation (j5.13p is the 

non-conforming error and that it is bounded by C*!! Vn||(j;^2(Q))d sup^^gy ""llcil' 

an appropriate constant C). 

l|Vn-C||(L2(Q))d < C(inf ||Vn-C||(L2(t^))d + || Vu||(i2(Q))d sup ^l^^J^j^!!^^-^) (5.17) 

Combining (|5.16|) with (|5.17|) . we conclude using theorem l2.2l and the Poincare inequality. 

□ 

Let us now show how theorem 15.11 can be combined with the new class of inequalities 
obtained in sub-section 14. 1 1 to obtain homogenization results for arbitrarily rough coeffi- 
cients a. Let M be a uniformly elliptic d x d matrix (observe that uniform ellipticity of 
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M implies its invertibility) and V be a finite dimensional linear subspace of (L^(O))'^. 
Define 

V := {MC : C' G V'} (5.18) 
Assume furthermore that for all w E Hq H H'^{Q.), 

^inf ^ ||Vu; - C'II(L2{Q)) < C/i||Au;||i2(n) (5.19) 

where /i is a small parameter (the resolution of the tessellation associated to V' for 
instance). We remark here that V' can be viewed as the coarse scale h approximation 
space (see example below). The fine scale information from coefficients a(x) is contained 
in the elements of the matrix M. This is illustrated in the example below where M = VF 
for harmonic coordinates F. Therefore, the matrix M is determined by d harmonic 
coordinates that are analogues of d cell problems in periodic homogenization, and we 
call space V the "minimal pre-computation space" since it requires minimal (namely d) 
pre-computation of fine scales. 

Then, we have the following theorem: 

Theorem 5.2. Approximation by "minimal pre-computation space" // 

• a ■ M is divergence free (as defined in sub-section\4.1\). 



• The symmetric part of a-M satisfies the Cordes condition ( |^.^7 ) or the symmetric 
part of a ■ M satisfies one of the inequalities of theorem 

• The non- conforming error satisfies fCy < C/i" for some constant C > and 
a G (0,1], 

then 

\\^u-CvkL^n)Y<C\\f\\mn)h'' (5.20) 
where u is the solution of and Qy the solution of \5.^) . 

Remark 5.4. The error estimate is given in the norm of Vn — Cv because we wish to 
give a strong error estimate and if is not a gradient then the norm of {a{Vu—(!^\;))pot 
is not equivalent to the norm Vu — (^y 

Remark 5.5. It is, in fact, sufficient that the symmetric part of a • M satisfies one of the 
inequalities of theorem 14.11 instead of 14.21 for the validity of Theorem 15.21 For the sake 
of clarity, we have used inequalities of theorem] 



Proof. The proof is a direct consequence of theorem I5.lt we simply need to bound Dy. 
Since div(aV) = div(a • MV'), it follows from equation 15.41 that 



\{a-M{Vu' -C))pot \\ (i2(n))rf 



Vv < sup inf — — "^"^ ^"^^ (5.21) 

- ^ A/^i,/ I At(;'ll-"— ^ ^ 
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where n' in (|5.4|) is defined as the (weak) solution of div(a-MVii') = Aw' with Dirichlet 
boundary condition on d^l. Now, if symmetric part of a-M satisfies the Cordes condition 
(j4.27p or the symmetric part of a • M satisfies one of the inequaHties of theorem 14.21 then 
1 1 1 1/^2 ^ CII^^^'IIl^ and we conclude using the approximation property ()5.19p . 

□ 

An example of V can be found in the discontinuous Galerkin method introduced 
in subsection 1.3 of |53j . This method is also a generalization of the method II of 
[8] to non-laminar media. In that method, we pre-compute the harmonic coordinates 
associated with (|l.ip -i.e.. the d-dimensional vector F{x) := . . . where Fj 

is a solution of 

div aVFi = in Q , , 

(5.22) 

Fi[x) = Xi on do,. 

Introducing Q^, a regular tessellation of of resolution h, the elements of V are defined 
as VF{VcF)~^Vip, where 93 is a piecewise linear function on with Dirichlet boundary 
condition on il.h and VcF is the gradient of the linear interpolation of F over Qh- In 
that example a ■ VF is divergence-free and VF plays the role of M. The non-conforming 
error is controlled by the aspect ratios of the images of the triangles of O/i by F. In [53], 
the estimate (|5.20p is obtained using F as a global change of coordinates that has no 
clear equivalent for tensorial equations, whereas the proof based on the flux- norm can 
be extended to tensorial equations. 



5.2 Tensorial equations. 

The generalization of the results of this section to elasticity equations doesn't pose any 
difficulty. This generalization is simply based on theorem 12.41 and the new class of 
inequalities introduced in subsection 14.21 An example of numerical scheme can be found 
in [38] for (non- linear) elasto-dynamics with rough elasticity coefficients. With elasticity 
equations harmonic coordinates are replaced by harmonic displacements, i.e. solutions 
of 

f-div(C7(x)VF*^0 = xefl 

and strains £{u) are approximated by a finite dimension linear space V with elements of 
the form e{F) : {ecF)~^{e{ip)) where the ip are piecewise linear displacements on ^1^, EcF 
is the strain of the linear interpolation of F over and £{F) denotes the d x d x d x d 
tensor with entries 

a^Ff + diFj'^ 

e{F).,,k,i ■■= ' ^ ■ (5-24) 

Here, C : £{F) is divergence-free and plays the role of M; furthermore, the regularization 
property observed in the scalar case [53] is also observed in the tensorial case by taking 
the product {£{F))~^£{u) (figured]). 



32 




(c) 

Figure 1: Computation by Lei Zhang. The elasticity stiffness is obtained by choosing its 
coefficients to be random and oscillating over many overlapping scales. Figure (a) and (b) 
show wild oscillations of one of the components of the strain tensor Vu + 'Vu'^ (u solves 
(fO]) ) and one of the components of {VF + VF^)-^ (F = {F^i} is defined by ([OH]) ). 
Figure (c) illustrates one of the components of the product {SI F + V F'^)~^{SIu + Vu'^), 
which is smooth if compared to (a) and (b). There is no smoothing near the boundary 
due to sharp corners. 

6 Relations with homogenization theory and other works. 

We first show how our approach is related to homogenization theory. To this end, we 

• Describe the notion of a thin subspace, which is pivotal in our work and show that 
this notion was implicitly present in classical periodic homogenization. 

• Show the analogy between the basis functions in Theorems 13.11 and 13. 3^ harmonic 
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coordinates (j4.35|) . and solutions of the cell problems in periodic and random ho- 
mogenization. 

• Explain relations between our work and general abstract operator liomogenization 
approaches. 

The thin subspace notion. A key ingredient of the proofs of the main approxi- 
mation Theorems 13.11 and 13.31 is the transfer property introduced in theorems 12.21 and 
12. 3[ Roughly speaking, it shows how a standard ("easy") error estimate in the case of 
smooth coefficients provides an error estimate in the case of arbitrarily rough ("bad") 
coefficients due to an appropriate choice of the finite-dimensional approximation space. 

The transfer property in turn is based on the notion of a "thin" subspace whose 
essence can be explained as follows. Let us consider the scalar divergence form elliptic 
problem First, observe that as / spans H~^{Q), u spans //(^(Jl), i.e. the operator 

:= (— divaV)"^ defines a bijection from H~^{Q) onto Hq{Q). Next, observe that 
as / spans L^(r2), u spans a subspace V of Hq{Q,), i.e. := (— divaV)^^ defines a 
bijection from L^(r2) onto V. How "thin" is that space compared to Hq{Q,)7 If a = Id, 
then V = V := Hq{Q) n H'^{Q), i.e. a much "thinner" space than Hq{^1,), namely V is 
"as thin as i?^" . 

The proofs of Theorems l3.1l and l3.3l also use the transfer property for finite-dimensional 
approximation spaces that depend of a small parameter h. Namely, there exists a finite 
0{\Q\/h'^) dimensional subspace of i?g(0) such that all elements of y' are in i7^-norm 
distance at most h from (an example of spaces and V are the spaces Cq (used in 
subsection 13. l.ip and H'^{Q) respectively). 

Section [3] shows that when the entries of a are only assumed to be bounded, the 
solution space V is isomorphic to V' = Hq{Q) n H'^{Q), that is for arbitrarily rough 
coefficients the approximation space V is still "as thin as" (isomorphic to H^). 
Moreover, the transfer property introduced in theorem 12.21 allows us to explicitly con- 
struct a finite-dimensional space V^, isomorphic to V^, such that all elements of V are 
in norm distance at most h from Vh- 

We next show that the thin subspace notion is implicitly present in classical homog- 
enization problem when a is periodic, with period e, i.e. when equation (jl.ip is of the 
form 

- diY{a{x/e)Vu^x)) = f{x), inQcR'^ (6.1) 

From the two-scale asymptotic expansion ansatz justified in periodic homogenization 
(e.g., [13], [12], [23], [3]), we know that can be approximated in norm by 

(modulo boundary correctors which we do not discuss here for the sake of simplicity of 
presentation) 

where u is the solution of the homogenized problem 

- div(aVu) = /(x) (6.3) 
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with constant homogenized (effective) coefficient a. Here, the exact solution has both 
fine, 0(e), and coarse, 0(1), variations (oscillations), while the homogenized solutions 
u{x) has only coarse scale variations. In (j6.3p . the periodic functions Xk a^^e solutions of 
the cell-problems 

div(a(y)(efc+Vxfc(y))) =0 (6.4) 

defined on the torus of dimension d. The second term in the right hand side of (|6.3p is 
known as a corrector; it has both fine and coarse scales, but the fine scales (e- oscillations) 
enter in a controlled way via d solutions of the cell problems that do not depend on f{x) 
and the domain (so that Xk are completely determined by the microstructure a{x)). 

Furthermore, the two-scale convergence approach jS^ [2] provides a simple and el- 
egant description of the approximation for the gradients Vn^. Namely, for every suffi- 
ciently smooth 4){x, y) which is also periodic in e, ^ u weakly in i?^(i7) and 

/" 0(x, -)Vn, ^ /" <i){x,y){Id + Vx.{y))^u{x) (6.5) 

where T"' is the torus of dimension d. Thus Vu^ can be approximated in the sense of 
()6.5p by functions of the form [Id + Vx.{^/^))^'^{^) 

The latter observation combined with (j6.2p shows that classical homogenization re- 
sults can be viewed as follows. The solution space V for the problem ()6.ip can be 
approximated in norm by a ("thin") subspace of H^{Q) parameterized by solutions 
of (|6.3p which are in n Hq{Q,). Moreover, homogenization theory shows us how to 
construct an approximation of the space V. Indeed, (|6.2p . (|6.5p show that this approx- 
imation space is determined by u{x) and solutions of the cell problems ()6.4p over one 
period. 

Cell problems. This periodic homogenization scheme was generalized to stationary 
ergodic coefficients a{x/e,uj), with ui in some probability space. Here, there are also 
analogs of the cell problems that require the solution of d different boundary value 
problems for the PDE div {a{x,uj)Vui) = 0, i = 1, . . . , d in a cube of size i? — t- oo for a 
typical realization oj |371I56|. The solutions Ui must be pre-computed in order to obtain a 
homogenized PDE and an approximate solution just like in the periodic problem, which 
is why they still are called cell problems even though there is no actual periodicity 
cell in the microstructure. Here, the coefficients in the cell problems have both fine 
and coarse scales, as in the original PDE div{a{x/e,uj)'Vu^) = f. The advantage of 
solving cell problems numerically (in both the periodic and random case) comes when, 
e.g., we need to solve for many different / or for a corresponding evolution problem 
= div(a(x/e, w)Vti^) — / when updating in fine time scales. 

Note that the major difficulty in advancing from periodic to random homogenization 
was to understand what is the proper analog of the periodic cell problem. In this work, we 
ask a similar question- what are the analogs of cell problems for most general arbitrarily 
rough coefficients? For the problem (jl.ip . we provide two answers. 

First, in Theorem 13. 3t we introduce functions 6k{x) that are analogs of the cell prob- 
lems since they are determined by the coefficients a{x) but do not depend on f{x). 
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Note, that these "generalized cell problems" must depend on the domain since the 
coefficients are no longer translationally invariant (as in periodic and random stationary 
case). This dependence enters via in (|3.8p using the transfer property. The approx- 
imation space in Theorem 13.31 is Qh defined in ()3.9p - ()3.f0p . Similarly, in Theorem 13. 1| 
we introduce functions ^ki^) that solve (j3.ip with localized right hand side that are also 
analogs of the cell problems. 

Secondly, the harmonic coordinates (j4.35p provide yet another analog of cell problems 
in classical homogenization. Their advantage is obvious- there is only d of them, whereas 
in Theorems 13.31 and 13.11 the number of cell problems (number of elements in the basis 
of Vh) is N{h). In fact, since in the simplest case of periodic homogenization d cell 
problems ()6.4p must be used, one should expect that for the more general coefficients 
a{x), d would be the minimal number of cell problems. On the other hand, the finite 
dimensional approximation based on harmonic coordinates, in general, is not direct 
(involves the non-conforming error) as one can be seen from Theorem 15.21 and example 
right after this Theorem. 

Harmonic coordinates play an important role in various homogenization approaches, 
both theoretical and numerical, which is why we present here a short account of their 
development. Recall that harmonic coordinates were introduced in [H] in the context of 
random homogenization. Next, harmonic coordinates have been used in one dimensional 
and quasi-one dimensional divergence form elliptic problems [101 [8] , allowing for efficient 
finite dimensional approximations. 

The idea of using particular solutions in numerical homogenization to approximate 
the solution space of (jl.ip have been first proposed in reservoir modeling in the 1980s 
|18j . |63j (in which a global scale- up method was introduced based on generic flow solu- 
tions i.e., flows calculated from generic boundary conditions). Its rigorous mathematical 
analysis was done only recently [53]. In [53], it was shown that if a{x) is not periodic 
but satisfies the Cordes conditions (a restriction on anisotropy for d > 3, no restriction 
for d = 2 and convex domains $7), then the ("thin") approximation space V can be 
constructed from any set of d "linearly independent" solutions of (jl.ip (harmonic coor- 
dinates F, for instance by observing that uo spans n Hq{Q,) as / spans L^(J7)). 
In the present work (section [5|) for elasticity problems, we show that d{d + l)/2 "linearly 
independent" solutions are required (equation (|5.23p ). 

In [3], the solution space V is approximated by composing splines with local harmonic 
coordinates (leading to higher accuracy) , and a proof of convergence is given for periodic 
media. Harmonic coordinates have been motivated and linked to periodic homogeniza- 
tion in |3] by observing that equation ()6.2p can in fact be seen as a Taylor expansion 
of u{x + x-i^)) where x + xX^) is harmonic, i.e., satisfies ()6.4p . It is also observed in 
[3] that replacing x + x.{^) by global harmonic coordinates F{x) automatically enforces 
Dirichlet boundary conditions on -u^. 

More recently, in |29l [281 dS] > idea of a global change of coordinates analogous to 
harmonic coordinates was implemented numerically in order to up-scale porous media 
flows. We refer, in particular, to a recent review article [18] for an overview of some main 
challenges in reservoir modeling and a description of global scale-up strategies based on 
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generic flows. 

Abstract operator homogenization approaches and their relation to our work. 

Recall that the theory of homogenization in its most general formulation is based on ab- 
stract operator convergence, -i.e., G-convergence for symmetric operators, i^-convergence 
for non-symmetric operators and F-convergence for variational problems. We refer to 
the work of De Giorgi, Spagnolo, Murat, Tartar, Pankov and many others [HI [311 [32l 
[591 l58l H71 [T7]). H, G and F-convergence allows one to obtain the convergence of a se- 
quence of operators parameterized by e under very weak assumptions on the coefficients. 
The concepts of "thin" space and generalized cell problems are implicitly present in 
this most general form of homogenization theory through the introduction of oscillating 
test functions in ff-convergence |48] (see also related work on G-convergence |581 131j). 
Furthermore, the so called multiscale finite element method |351 [6l] can be seen as a 
numerical generalization of this idea of oscillating test functions with the purpose of con- 
structing a numerical (finite dimensional) approximation of the "thin" space of solutions 
V . We refer to [33] for convergence results on the multiscale finite element method in 
the framework of G and F-convergence. 

Observe that in most engineering problems, one has to deal with a given medium and 
not with a family of media. In particular, in those problems, it is not possible to find a 
small parameter e intrinsic to the medium with respect to which one could perform an 
asymptotic analysis. Indeed, given a medium that is not periodic or stationary ergodic, it 
is not clear how to define a family of operators A^. Moreover, the definition of oscillating 
test functions involves the limiting (homogenized) operator A. While this works well for 
the proof of the abstract convergence results, in practice only the coefficients A are 
known (computing A may not be possible), and our approach allows one to construct 
the approximate (upscaled) solution from the given coefficients without constructing 
A. Hence, the main difference between H or G convergence and present work is that 
instead of characterizing the limit of an e-family of boundary value problems we are 
approximating the solution to a given problem with a finite-dimensional operator with 
explicit error estimates (i.e. constructing an explicit finite dimensional approximation 
of the "thin" solution space V). 

We refer to [25] for for an explicit construction of A with rough coefficients a in 
two dimensions. In particular, it is shown in [25] that conductivity coefficients a are 
in one-to-one correspondence with convex functions s(x) over the domain Q and that 
homogenization of o is equivalent to the linear interpolation over triangulations of 
re-expressed using convex functions. 

The thin subspace idea introduced in this section can be used to develop coarse 
graining numerical schemes through an energy matching principle. We refer to [S5] for 
elasticity equations and to |65j for atomistic to continuum models (with non-crystalline 
structures). 

Other related works. By now, the field of asymptotic and numerical homogenization 
with non periodic coefficients has become large enough that it is not possible to cite all 
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contributors. Therefore, we will restrict our attention to works directly related to our 
work. 

- In [26\ I30j . the structure of the medium is numerically decomposed into a micro- 
scale and a macro-scale (meso-scale) and solutions of cell problems are computed on the 
micro-scale, providing local homogenized matrices that are transferred (up-scaled) to the 
macro-scale grid. This procedure allows one to obtain rigorous homogenization results 
with controlled error estimates for non periodic media of the form a{x, |) (where a(x, y) 
is assumed to be smooth in x and periodic or ergodic with specific mixing properties 
in y). Moreover, it is shown that the numerical algorithms associated with HMM and 
MsFEM can be implemented for a class of coefficients that is much broader than a{x, ^). 
We refer to jSSj for convergence results on the Heterogeneous Multiscale Method in the 
framework of G and T-convergence. 

- More recent work includes an adaptive projection based method [51], which is 
consistent with homogenization when there is scale separation, leading to adaptive al- 
gorithms for solving problems with no clear scale separation; fast and sparse chaos 
approximations of elliptic problems with stochastic coefficients [611 ; finite difference 
approximations of fully nonlinear, uniformly elliptic PDEs with Lipschitz continuous 
viscosity solutions [21] and operator splitting methods [71 [6] . 

- We refer the reader to [60] and the references therein for a series of computational 
papers on cost versus accuracy capabilities for the generalized FEM. 

- We refer to [W\ [TS] (and references therein) for most recent results on homoge- 
nization of scalar divergence-form elliptic operators with stochastic coefficients. Here 
the stochastic coefficients a{x/e,u}) are obtained from stochastic deformations (using 
random diffeomorphisms) of the periodic and stationary ergodic setting. 

- We refer the reader to [22] , [27] and |14j for recent results on adaptive finite element 
methods for high contrast media. Observe that in [22], contrast independent error 
estimates are obtained for a domain with high contrast inclusions by dividing the energy 
norm by the minimal value of a over the domain il. The strategy of [H] is to first prove 
a priori and a posteriori estimates that are contrast independent and then construct a 
finite element mesh adaptively such that the error is the smallest possible for a fixed 
number of degrees of freedom. In [14] . the energy norm of the error is bounded by terms 
that are appropriately and explicitly weighted by a to obtain error constants independent 
of the variability of a. 

7 Conclusions 

In this paper, we have been primarily concerned with the following theoretical results: 

• the flux norm introduced in section [2] and its transfer property (theorems 12.21 and 

1231) 

• the resulting theoretical Galerkin method with contrast independent accuracy for 
linear PDEs with arbitrarily rough coefficients (theorems 13. H 13.31 13.41 13.30p 
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• a new class of inequalities (section H]) and the introduction of a Cordes-type con- 
dition for non-divergence tensorial elliptic equations (theorem 14. 6p . 

The development of numerical techniques based on these results is in progress and will 
be addressed elsewhere. 



8 Appendix 

8.1 Extension to non-zero boundary conditions. 

The analysis performed in section [3] and the previous one can naturally be extended to 
other types of boundary conditions (Neumann or Dirichlet). To support our claim, we 
will provide here this extension in the scalar case with non-zero Neumann boundary con- 
ditions. By linearity, solutions of (jl.ip with nonzero boundary conditions can be written 
as the sum of the solution with zero-boundary condition and — div (^a{x)Vu{x)^ = f{x) 

and the solution with non-zero boundary condition and — div ^a(x)Vu(x)^ = 0. Hence, 
we will restrict our analysis to solutions of 



We assume d > 2. Write L'^otni^) the closure of {Vv : v G H^{^)} in {L^{n)f. Note 
that the difference with Lp^f{Q) lies in replacement of G //g(rj) hy v € H^{i}) in 
the definition of Lp^^j^{Q). For ^ in (L^(r2))'^, write S,potn its orthogonal projection on 
L%,^{n). For V G H\Q), write 

ll^'lliflux,!! := / {0'^v)l^tn (8-2) 

It follows from the following lemma that the ||u||a_flux,n-noi'™ is equivalent to the H^- 
norm. 

Lemma 8.1. For v G H^{Q) 

Amin(a)||V7;||(i2(f^))d < ||u||a-flux,n < Vax (a) || Vu || (^2 (8.3) 



Proof. The proof of the right hand side of (|8.3p is straightforward. The proof of the left 
hand side follows from 

j (V?;)'^aV?; < ||V?;||(i2(!^))d||i;||a_fiux,n (8.4) 

□ 
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Write A the Dirichlet to Neumann map mapping v\qq onto n • Vv on dil. where v is 
the solution of Av = in (and n, the exterior normal to the boundary dQ). Write 

the orthonormal eigenvectors of A and the associated increasing and positive 
eigenvalues. Let Vn be the finite dimensional subspace of H^{Q,) formed by the linear 
span of vi, . . . ,Vn, where Vk is the solution of (jS.ip with h = '^f^. 

Theorem 8.1. For h G L^{Q), let u be the solution of Ii8. 1]) . Then, 

sup mf — —r. = r — 0.5 

Furthermore, the space Vn leads to the smallest possible constant in the right hand side 
of 118. 5|) among all subspaces of H^{^) with n elements. 

Remark 8.1. We also refer to [2] for the related introduction of the penetration function, 
which measures the effect of the boundary data on the energy of solutions of (|8.ip and 
is used to assess the accuracy of global-local approaches for recovering local solution 
features from coarse grained solutions such as those delivered by homogenization theory. 

Proof. For h G L^{dQ), let w'^ be the solution of 
(Aw^ix) = xGQ; 

\n-Vw'^'{x) = h{x) xedn; h£L^{dn), 

Let Wn be a finite dimension subspace of H^{Q,) of elements v satisfying div(aVf) = 
in $7. To prove Theorem 18.11 we will first prove lemmas 18.21 and 18.31 given below. 

Lemma 8.2. For h G L'^{dQ,), let u be the solution of \8. 1\) and be the solution of 
[Km . We have, 

. r \\U - V\\a-fln^,n . , II (Vw^ - aVujpotn || (^2 (n))d 

sup mf — — = sup mt — (8.7) 

Proof. The proof follows by observing that for all ip G H^{il) 

/ VifiVw^' - aVu) = / tpn- (yw^' -aVu) = (8.8) 
Jn Jan 

□ 

For V G H^{^), write h^ := n.aVv defined on dQ 
Lemma 8.3. For h G L^{dn) and v G H^{^), 

\\{Vw'' - aVv)potn\\lL2in)y = [ {h - h^)k-\h - (8.9) 
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Proof. The proof is obtained by observing that 



h ^ „ LVifiVw'^ - aVv) 
\\[Vw ■ - aVv)potn\\(L^n)y = sup -— 

— sup 



/^VyVit;(^~^") (8.10) 

sup 1 1 j"j 

{h - h'')A-\h - h") 

an 

□ 

Write Zn the subspace of H~2(^dQ) induced by elements n ■ aVv defined on dQ for 
V € Wn- Write /i^ the coefficients of h in the basis i-e- ^fc := JoQ^^ki It follows 
from lemmas 18.21 and 18.31 that 



. r Ik - ■^Ila-fiux,n . , (YX=1 Ti^^^k - Kf' \\ 

sup mf — — = sup mf =^55 — -3 (o-H) 

It follows that an space Wn with optimal approximation constant is obtained when Z„ 
is the linear span of (^1, . . . , ^„,). In that case 

■ f 11^ ~ ^l|a-flux,n 1 /'0 10^ 

sup mf — — = Y (8.12) 

/iei:2(aQ)^ew„ \\n\\L-2{dQ) (A„+i)2 

We deduce equation ()8.5p using the following Weyl's asymptotic for A (as a first 
order pseudo-differential operator, see for instance [30], [l3] and 



\dVl\ 

This concludes the proof of Theorem 18.11 



(Afc)^^4vr( ^^^ + 'pS --\ (8.13) 



□ 



8.2 Proof of lemma 14.11 

Let u be the solution of Cu = f with Dirichlet boundary condition (assume that it 
exists). Since a > 0, the solvability of ()4.3ip is equivalent to finding u G H Hq{Q,) 
such that 

Au = af + Au- aCu (8.14) 
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Consider the mapping T : n Hq{^) — )• i/^ Pi Hq{^) defined by u = Tw where v 
be the unique solution of the Dirichlet problem for Poisson equation 

^v = af + ^w-aCw (8.15) 

Let us now show that for < 1, 7" is a contraction. 



||Tw;i - Tu;2||^2n^i(^) = - V2\\H2r^Hl{n) (8.16) 

Using the convexity of one obtains the following classical inequality satisfied by the 
Laplace operator (see lemma 1.2.2 of j45j ) 

\\vi - V2\\mr)Hl{n) < W^i'^i - V2)\\L^{n) (8.17) 

Hence, 

||ru;i - Tw2\\\2^Hl{Q.) ~ ^2) - aC{wi - W2)\\\2^a-^ 

2 (8.18) 



-| X] {^^i ~ aaij)didj{wi - W2) 



L2(n) 



Using Cauchy-Schwarz inequality we obtain that 

d 



\Twi - Tw2\\l2nH^^^n) - iYl ^^^^ " 



; 5^ (9.9,(^-4))') 



.19) 



Hence observing that 

d 

esssupQ( ^ {5ij - aatjf) = Pa (8.20) 

we obtain that 

\\Twi - Tw2\\\2^Hl(n) ^ esssup^gf2/3a(x)||u;i - «^2||^2n/^i(n) (8-21) 

It follows that if /3c < 1, then T is a contraction and we obtain the existence and solution 
of ()4.3ip through the fixed point theorem. Moreover, 

II Au|U2(f^) < ||a/||L2(o) + /3| II An||i2(t^) (8.22) 
which concludes the proof. 
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